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Abstract 

Li-Vogelius and Li-Nirenberg gave a gradient estimate for solutions of 
strongly elliptic equations and systems of divergence forms with piecewise 
smooth coefficients, respectively. The discontinuities of the coefficients are 
assumed to be given by manifolds of codimension 1 , which we called them 
manifolds of discontinuities. Their gradient estimate is independent of the 
distances between manifolds of discontinuities. In this paper, we gave a 
parabolic version of their results. That is, we gave a gradient estimate for 
parabolic equations of divergence forms with piecewise smooth coefficients. 
The coefficients are assumed to be independent of time and their disconti- 
nuities are likewise the previous elliptic equations. As an application of this 
estimate, we also gave a pointwise gradient estimate for the fundamental so- 
lution of a parabolic operator with piecewise smooth coefficients. The both 
gradient estimates are independent of the distances between manifolds of 
discontinuities. 



1 Introduction. 

For strongly elliptic, second order scalar equations with real coefficients, it is 
well-known that their solutions have the Holder continuity even in the case that 
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the coefficients are only bounded measurable functions. However, the solutions do 
not have the Lipschitz continuity in general. For example, Piccinini-Spagnolo IfTSl 
p. 396, Example 1] and Meyers fil2L p. 204] gave the following example: 

Example 1.1. (Ea, [USD Let Bi := {x G M" : |x| < 1} and each G L°°{Bi) 
be defined as 

_ Mxl + xl _x\ + Mxl _ _ (M - l)xiX2 

C^ll — i — 5 ^22 — j — 77, , ai2 — 021 — j — 77, 

with a constant M > 1. Then, if we define u as 

(1.1) 



U[X] = \x\ 



\x\ 



it is easy to see that the Holder exponent of u is at least less than or equal to 1/ v M 
(indeed, for a; = (xi, 0) we have \u(x) — m(0)| = \x\^/^ . Hence we have 

' ^ ^ ^ '-\x\ ^ ^ +00 as a; 



\x\ 



(l/v^)+£ 



for any e > 0.) and u satisfies the strongly elliptic scalar equation with real 
coefficients 



The same thing can be said also to the parabolic equation 

du d f du\ ^ 

dt ^ dxi \ dxi J 

i,j=l ■' 

because u given by (11.11) satisfies this equation. 

This example shows that we cannot expect gradient estimates of solutions to 
equations (11.21) and (|1.3I) in the case ajj G (i^i), but we may have the estimates 
in the case of piecewise C^' (see (|1.5I) below) coefficients. 

The fact that the gradient estimate of solutions is independent of the distances 
between manifolds of discontinuities was first observed by Babu ska- Anders son- 
Smith-Levin [2] numerically for certain homogeneous isotropic linear systems of 
elasticity, that is |Vm| is bounded independently of the distances between mani- 
folds of discontinuities. They considered that this numerical property of solutions 
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is mathematically true. This is the so-called Babuska's conjecture. Recently, [fTTTl 
and fTO'l gave mathematical proofs for this conjecture. In elasticity, a small static 
deformation of an elastic medium with inclusions can be described by an elliptic 
system of divergence form with piecewise smooth coefficients. The discontinu- 
ities of coefficients form the boundaries of inclusions. Similar physical interpre- 
tation is also possible for heat conductors. Our main theorem [1.51 given below 
ensures that this property also holds for parabolic equations of the form (11.31 ). The 
details of result given in [1 IJ and filOJ for scalar equations will be given below as 
Theorem 1 1.2 1 

In order to state our main theorem, we begin with introducing several no- 
tations which will be used throughout this paper. Let D C R" be a bounded 
domain with a C^ " boundary for some < a < 1, which means that the do- 
main D contains L disjoint subdomains Di, . . . , Dl with C^ " boundaries, i.e. 
D = (IJm=i ^rn) \ dD, and we also assume that Dm C D for 1 < m < L — 1. 
Physically, D is a material and Dm (1 < m < L — 1) are considered as in- 
clusions in D. We define the C^'" norm (resp. C^ " seminorm) of C^'" domain 
Dm in the same way as in [[TOll . that is, as the largest positive number a such 
that in the a-neighborhood of every point of dDm, identified as after a possi- 
ble translation and rotation of the coordinates so that x„ = is the tangent to 
dDm at 0, dDm is given by the graph of a C^'" function ipm, defined in \x'\ < 2a 
(x' = (xi, . . . , Xn-i)), the 2a-neighborhood of in the tangent plane, and it satis- 
fies the estimate ||^/'m||ci>-(|x'|<2a) < l/a (resp. [^m]c^,-'{\x'\<2a) < whcrc 



Further, let (ajj) be a symmetric, positive definite matrix-valued function de- 
fined on D satisfying 




IIV^IIci>'»(|x'|<2a) ||^||cl(|2^'l<2a) + bP]c^''^ {\x'\<2a) ■ 



n 




(1.4) 



Here each aij is piecewise in D, < /i < 1, that is 



ajj(x) = ai^ \x) for a; G -Dm, I < m < L 



(1.5) 



witha^^ eC^(Dm). 
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As we have already mentioned above, we will discuss in this paper a gradient 
estimate for solutions to parabolic equations with piecewise smooth coefficients. 
Our result is a parabolic version for the results of Li-Vogelius [1 1 J and the scalar 
equations version of Li-Nirenberg [10 J . They showed that solutions u E H^{D) 
to the elliptic equation 

jj=i ^ ' 1=1 

where h E L°°{D) and each gi is defined in D such that giln^ (1 < < -^^) have 
continuous extensions E C^{Drn), < < 1 up to dDm have global W^'°° and 
piecewise estimates (see (11.71) below). These estimates are independent of 
the distances between inclusions when a material has inclusions. 

We first give the result of Li-Nirenberg [|TOl for scalar equations. 

Theorem 1.2 ( [[TOl Theorem 1.1]). For any e > 0, there exists a constant Cjj > 
such that for any a' satisfying 

I f 

< a < mm < /i 



'2(a + ll 



we have 

L / L n 

m=l \ m=l 1=1 

(1.7) 

where we denote 

De:= {xE D : dist(x, dD) > e} 

and a positive constant depends only on n, L, fj,, a, e, A, A, ||aij||ca'(z);;;) (^^^d 
the C^'"' norms of Dm- 

Remark 1.3. The constant Q > is independent of the distances between in- 
clusions Djn. Therefore, the estimate (11.71 ) holds even in the case that some of 
inclusions touch another inclusions as in Figure [U 

Now, we consider the parabolic equation 

du ^ d ( du\ ^dfi 

1,3=1 ■' J=l 
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Figure 1: The case that an inclusion touches another inclusion. (L = 7) 



where 



h G U{Q), ^ G U{Q) and /, = on x (0, T], 



with p > n + 2, K = p{n + 2)/{n + 2 + p), Q := D x (0,T] and /^""^ G 
L-(0,T;C'^(A^)). 

Now we define a weak solution to the equation (|1.8b . 

Definition 1.4. We call u G l^2^'°(Q) := L\0,T;H\D)) n C([0, T]; L^p)) a 
weak solution to the equation (11.81) when 

f^' f dip 
u{x,t') ip{x,t') dx — / / u{x,t) —{x,t) dx dt 



D 







dt 



+ 



JD 



dxdt+ [ [ V/i(x,t) |^(x,t)dxdt (1.9) 

Jo JD (J^i 



f{x,t)ip{x,t) 



for any G L^{0,T; H\D))nH\0,T; L^{D)) with ip{-,0) = and < t' < T. 
Our main result is as follows. 
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Theorem 1.5 (Main theorem). Any weak solutions u G V^2^'°(<5) to (11.81) have the 
following up to the inclusion boundary regularity estimate: For any £ > 0, there 
exists a constant Cu > such that for any a' satisfying 



< a' < mill < fi 



a 



2{a + l) J ' 



(1.10) 



we have 



^ sup ||M(-,t)|lciV(s;;;nD,) < C'fl (||m||l2(q) + F, + 



m=l 



£'^<t<T 



where 



df 



dt 



F„ :- ^ ||/i|Up(Q) + 



1=1 



dfi 



dt 



LHQ) 



^f^ 



dt 



LP{Q) 



m=l 



0<t<T 



0°^ (Dm) 



and C'^ depends only on n, L, /i, a, e, A, A,p, ||oij||(7"'(z);;r) ^"'^ C^'"' norms of 

Dm. 



Remark 1.6. (i) Again, the constant > is independent of the distances be- 
tween inclusions D^- Then Theorem II 51 holds even in the case that an inclusion 
touches another inclusion as Figure [1] 
(ii) It is easy to obtain 



F, < C* 





df 




y\f\\L^iQ) + 


dt 





n / L 



However, a constant C* > depends on T and D, unfortunately. 







dt 
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For heat conductive materials with inclusions, (11.81 ) describes the tempera- 
ture distribution in the materials. When these inclusions are unknown and need 
to be identified, thermography is one of non-destructive testing which identifies 
these inclusions. The measurement for the thermography could be temperature 
distribution at the boundary generated by injecting heat flux at the boundary. The 
mathematical analysis for this thermography has not yet been developed so far. 
However, if we have enough measurements, the so called dynamical probe method 
(lEII) can give a mathematically rigorous way to identify these inclusions. In the 
proof of justifying this method, the gradient estimate of the fundamental solution 
of parabolic equation with non-smooth coefficient is one of the essential ingredi- 
ents. 

The dynamical probe method has been developed only for the case that the 
inclusions do not touch another inclusions. So, it is natural to consider the case 
when some of them touch. For the first task to handle this case, we need to have 
the gradient estimate of the fundamental solution. Our main result has given the 
answer to this. Similar situation can be considered for stationary thermography 
and non-destructive testing using acoustic waves. For example, lfT4l and lfT6ll 
effectively used a result of Li-Vogelius [illl to give a procedure of reconstructing 
inclusions by enclosure method (see for example). What is interested about 
their arguments is that, by adding further arguments, we can even reconstruct the 
inclusions in the case that they can touch another inclusions ( llT3l ). Therefore, we 
believe that our gradient estimates will be useful for inverse problems identifying 
unknown inclusions. 

The rest of this paper is organized as follows. In Section |2l we prove our 
main theorem, i.e. Theorem 11. 51 by applying Lemma IITI We prove Lemma \2A] 
in Section [31 In Section IH we consider a pointwise gradient estimate for the 
fundamental solution of parabolic operators with piecewise smooth coefficients 
by applying Theorem [T31 

2 Proof of main result. 

In this section, we prove our main theorem. We first state some estimates in 
Lemma [2T| which we need to prove our main theorem. We prove Lemma IXTI in 
Section |3l 

Lemma 2.1. Let (aij) be a matrix-valued function defined on D. Assume that 
{atj) is symmetric, positive definite, and satisfies the condition (11.41) . Let Q as 
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before and := x [e^, T]. Then for p > n + 2, a weak solution u G 
to (11.81 ) satisfies the following estimates: 



sup \\u{-,t)\\L2(D,) < C {\\u\\l^Q) + Fq) 



e'^<t<T 



du 



dt 



<C \\u\ 



+ F,) 



(2.1) 
(2.2) 
(2.3) 



where we set 

Fn := 



p(".+2) + 7 J\fi\\LP{Q), 



^1 '■= 11/11 max|2,P(" + ^)l + ll/i|Uf(Q) + 



5/. 



9t 



(2.4) 
(2.5) 



anJ C > depends only on n, A, A, p and e. 



Now we prove our main theorem by applying Lemma 12.11 This proof is in- 
spired by [jHl. 

ofTheorem \1.5\ Before going into the proof, we remark that a general constant 
C which we used below in our estimates depends only on n, A, A,p and Sj {j = 
1, 2, 3). To begin with the proof, let < < £2 < £^3- Then we have 



sup ||M(-,t)||L2(B^^) < C (||m||l2(q) +Fo) 



e^<t<T 



and 



du 



dt 



<C \\u 



+ F^) 



(2.6) 



(2.7) 



by (1211) and in Lemma O where Fq, Fi are defined by (l24l) and (1231) . On 
the other hand, Ut = du/dt satisfies the equation 



dut 
'dt 



d_ 

dxi 



aij{x) 



dut 
dxj 



df ^ d fdfi 



dt 



dxj \ dt 



i,j=l ' ^ -J ^ i=l 

by applying d/dt to (11.81) (also see Remark |Z2|) . Hence we have 



(2.8) 



by Lemma [2?n (12 .21 ). where we define 

df 



F' - 



dt 



p("+2) 



dt 



LPiQ) 



In particular, Ut(-,t) E L°°{Di;.J holds for a.e. t E {el,T]. Now we regard the 
equation (|1.8|) as the elliptic equation 



i,j=l ■' 



i=l 



(2.9) 



by fixing t E (e^, T]. We remark that dujdt - f E L°°(L>eJ. Then, for any a' 
with the condition (11.101) . we have the estimate 

L 

y^M^j- ^t)\\c^-"' (D^^nD..) 



m=l 



<cA \H-,t)\\L2^D.^) + 



du 



+ ||/(-,t)||L-{D,,) 
L n 

+ EEII/" 



(2.10) 



m=l j=l 



, ||"iillc"'{D„) 



by Theorem 11.21 where > depends only on n, L, fi, a, e, A, A 

and the C^'"' norms of Dm- Taking the supremum of the inequality (|2.10l) over 
{el, T] with respect to t, and using (|2.6I) . (|2.7I) and (|2.8I) . we have 

J2 sup ||M(-,t)|lci^"'(i5;;;nD,3) 



m=l ^2 



< CA sup ||M(-,t)i|L2(D,2) + 



,£|<i<T 



9t 



L n 

,^^P Il/i(-'^)llc;-'(i5;;;) 

m=l i=l £2<*<^ 



L n 



m=l i=l ^2 
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which is the estimate we want to obtain. 



□ 



Remark 2.2. Since we assume that u belongs only in V2 ' (Q) with respect to the 
regularity of a weak solution, one may think that we cannot apply d/dt directly. 
However, it is enough to consider the Steklov mean function and to make h tend 
to 0, where we define the Steklov mean function Vh of v by 

^ f-t+h 

Vh{x,t) = v{x,t) dr. 

Hereafter we omit the detail with respect to this remark although we often apply 
this argument. Also see flU III §2 p. 141] and (62) in (81 p. 152], for example. 



3 Some estimates. 

In this section, we prove Lemma 12.11 The estimates (12.11 ) and (12.21) are well- 
known, but we give these proofs in Appendix for readers' convenience. In order 
to show the estimate (12.31) . we prepare some necessary lemmas for its proof. 

Throughout this section, C > denotes a general constant depending only 
on n, A, A. Also, we assume that the coefficient (aij) is a matrix-value function 
defined on D, symmetric, positive definite, and satisfies the condition (11.41) . More- 
over, we set Qr := Br{xo) x (to — r^, to], and assume that Q2p C D x (0,T] with 
< p < 1. 

The following two lemmas are essentially shown in We give their proofs 
here for the sake of completeness. 

Lemma 3.1 ([8, Lemma 3]). Letl <r < 00 andl/r + l/r' = 1. Then a solution 
u to (11.81) satisfies the estimate 



Vm||l2(q^) < C 



1=1 



(3.1) 



Proof. Let ^ be a smooth cut-off function on Q2p satisfying ^ = 1 on Qp, C = 

on Q2p \ Qzp/2, < C < 1 on Q2p, and \dC/dt\ + | VCP < Cp-^ on Q^p. Let mq 
be the average value of u in Q2p: 

Uo '■= -rp\ — f // u{x,t)dxdt, 



\Q2p\ 



Q2p 
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where \Q2p\ denotes the measure of Q2p- Testing (11.81) by {u — Mo)C^ and inte- 
grating by parts (i.e. taking (p = {u ~ mo)C^ for (11-91 ). Also see Remark |Z2] ). we 
have 



1 
2 

+ 



B2p{xo) 



Edu du ^2 7 r 
dij — t; — C dxdt + 2 



Q2p 



[u — uoYC—dx dt 



y~] aij — {u - uo)C—dxdt 
Q2p^^ dxj 



' <92p 

Hence we have 



f{u — uq)C dx dt + 



2=1 



Q2p 



(/II 

ax,- ra 



dx dt. 



-! {{u-UofC^){x,to)dx + x[[ \Vu\^C^dxdt 



< 



1 



{{u - uoyC^){x,to) dx + 



B2p{xq) 

u — Uo)^(^dx dt — 2 



^2M 



Q2p 



+ 



dt 



f{u — Mo)C^ dx dt 



"^i^r^C dxdt 
"-ijTrri^ ~ Uo)C—dxdt 



' dxi 



Q2p 



i=l 



Q2p 



OXi OXi 



dx dt 



< II (u-uon 

'Q2p 



dt 



dxdt + Ei II \Vu\'^C'^ dx dt 

Q2p 



[[ \u-u,\^\VC?dxdt + \\ ff \fCYdxdt\ 

JjQ2p ^ \J JQ2p I 



c 
+ — 

^1 J JQ2p 



V 



1 

+ 2 



2/r' 



\{u-Ua)CY dxdt] \yu\X dxdt 



Q2p 



Q2p 



+ 



- + // Y,\u'edxdt+ [[ I 

^1 J J jQ2p J JQ2p 



\u - uq\^\WC\^ dxdt 
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We now take ei > small enough. Then, we have 



Vu\'^dxdt< II \Vu\'^Cdxdt 



<C 1 1 (u-uo) 

'Q2p 



c 



dt 



dx dt 



+ C 



2/r' 



\{u — 'Uo)Cr dxdt 



Q2p 



2/r 



< C 



+ C\II \fC\'dxdt\ +C ^l/ipC^rfxrft 

Q2p j J J Q2p i—\ 

2 n 

Q2p 



L^{Q2p) 

2=1 



because t) — Uq\ < oscqj^ u holds for any (x, t) G Q2p- This completes the 
proof. 

Lemma 3.2 ([^, Lemma 5]). A solution u to (|1.8I) satisfies the estimate 



□ 





< c 


at 





P l|VM||i2(Q2^) + \\f\\L\Q2p) 







dt 





(3.2) 



Proof. We first take the same smooth cut-off function C, as in the proof of Lemma lSTI 
Testing (11.81) by {du/dt)C'^ and integrating by parts (also see Remark IZ2] ). we have 



1 
2 

+ 



2 



9m (9m 

9xj dxj 



C {x, to) dx 



Q2p 



du 



dt 



du du ^dC, 
dxj dt dxi 



dx dt 



f^edxdt+Y: 

1=1 



dt 
+ 



B2p{xo) 



/i^C^ ) {x,to) dx 



Q2p 



dfi du ^2 o -C ^"^ /• , r, X ^"^ /■ "^C 



OT ax,- dxi dt dt dxj 



dx dt 
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due to 



dtdxi dxj 



^2J = 1 



du du 
dxi dxj 



i.j=\ 



du du ^dC, 
dxi dxi dt 



and 



Hence we have 



dtdxj 



d f „ du 



dt V dxj 



du d 
dxj dt 



^ f {\Vu\'e){x,to)dx+ [[ 

^ JB2p{xo) J JQ 

1 



< 





du 


/ 


dt 



C dx dt 



Edu du 2 \ / \ , 
I {x,to)dx + 



du 



Q2p 



^ du du dC ^ ^ 



Q2p 

du du ^ dC 



dx-i dt dxi 



dt 



dx dt 



dx dt 



,du 



Q2p 

+ 



1=1 



Jb2p(xo) V dXi J 



Q2p 



dfi du ^2 X- A^C , ox /■ 



OT dXj dXj dt dt dXi 



dx dt 
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<C [[ \Vu\\ 
J Jq^p 



dt 



dx dt + £2 



// 



du 



dt 



+ 



'Q2p 





du 


1 

IQ2p 


dt 



( dxdt 
2 

dx dt 



+ 



C 



II \f\'Cdxdt + e2 I {\Vu\^C){xMdx 

JjQin JBoJxn) 



+ 



£2 JJq^^ 

c 



£2 Jb2p(j;o) 



'B2p(a;o) 



.1=1 
I2/-2 



dt 



+ C If \VuW^dxdt + C If V 

JJQ2p J JQ2p j=l 

n 



Q"^ dx dt 



+ £2 
We remark that 



// 

J JQ2p 



du 



dt 



dt 



dx dt 



'dxdt+- II vi/invcr^xdi. 

£2 JJq,^ .^^ 



[ {M)\x,to)dx= [ r ^{{f,Cf){x,t)dtdx 

JB2p{xo) JBiJxn) Jtn-(2o')^ <^ 



< c 



B2p(xo) Jto-{2pY 

2 I a2 



Q2p 



dt 



+ 



Therefore, by taking £2 > small enough, we have 

(|V«|V) {x,t^)dx + 





du 






dxdt < 




dt 


J B2p(xo) 



dt 



du 



dx dt. 



<C [I |V«p 
J J 0.2. 



Q2p 



dt 



C dx dt 



2 / a2 



Q2p j=l 



dt 



+ |VCr ]dxdt + C II \ fix dxdt 

Q2p 



dt 



+ ivcr + 



dh 



dt 



dx dt 
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i=l 



dt 



LHQ2p) 



□ 



We obtain the estimate (12.31) from Lemmas lA. 5 1 (given in Appendix). 13.11 and 



4 A gradient estimate of the fundamental solution. 

In this section, we consider a gradient estimate of the fundametal solution of 
parabolic operators. We first state some facts. It is known that if coefficient (0,^) is 
a symmetric and positive definite matrix- value L°°(R") function satisfying (ll.4|) . 
then there exists a fundamental solution T(x, t; y, s) of the parabolic operator 




(4.1) 



with the estimate 

\T{x,t;y,s)\ < jj—^y^^^W ( — \_g j Xls,oo){t) (4.2) 

for all t, s G M, and a.e. x, y G M", where C*, c* > depend only on n, A, A (see 
ID or flU, for example). In particular, the constants C* and c* are independent of 
the distance between inclusions. If the coefficients (ajj) is not piecewise smooth 
but Holder continuous in the whole space M", then the pointwise gradient estimate 

/ c \x 1 ^ \ 

\V,T{x,t;y,s)\ < -^—-^-^exp xis,^){t) 

holds for t, s G M, a.e. x, ?/ G R" (see [9, Chapter IV §11-13], for example). 

Now, the aim of this section is to show the gradient estimate (14.81) in The- 
orem |43] even if the coefficients are piecewise in D. We assume that (aij) 
defined in D satisfies the conditions (|1.4I) and (11.51) . and extend it to the whole 
by defining (aij) = AI in W"- \ D, where / is the identity matrix. We remark that 



15 



this extension does not destroy the conditions (11.41 ) and (11.51) . Then there exists a 
fundamental solution T(x, t; y, s) of the parabolic operator (14. Il l with the estimate 
(14.21) as we stated above. 

To prove our gradient estimate of the fundamental solution, we apply the fol- 
lowing corollary from Theorem 11.51 

Corollary 4.1. Let < p < 1. Then a solution u to the parabolic equation 



du d f du\ ^ , s , 2 n . 

-Q^-l^Q^. y^'^d^. ) = Bp{xo) x{to-p , to] (4.3) 

i,j=l * ^ ^ 

has the estimate 

„ „ „ „ 

II^^IU°°(Bp/2(^o)x{to-{p/2)2,to])) < --;[J^\\u\\L^{B,ixo)x{to-p^,to]), (4-4) 

where (7^ > depends only on n, L, p, a, X, A, and ||c"'(I5;;r) '^^^ C^'"' 
norms of for some a' with (11.101) . 

Proof. It is enough to apply the scaling argument. To begin with, let py = x — xo, 

p^(s — 1) = t — to and 

u{y,s) := u{x,t) = u{py + Xo, p'^{s - 1) +to), (4.5) 
aij{y) := aij{x) = aij{py + xq), 

5m := I -(x - Xo) : X e D„ 
[P 

Then we have 

Therefore, by noting Remark 1421 we have 

II VM||ioo(Bj/2(o)x(3/4,l]) < C'^\\u\\l2(Bi{0)x{0,1)) 

by Theorem 11. 51 where C'^ depends only on n, L, p, a, A, A, ||ctij||c«'(i5;;;)' the 
(7^'"' seminorms of Dm- By this estimate and the definition (|4.51) . we obtain the 
estimate (l44l). □ 
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Remark 4.2. One may think that a constant depends also on p since || ctij || 

and the C^'"' norms of Dm depend on p. However, we can take C'^ independent of 
p by taking the following into consideration. 
First we consider 



:= sup \aij{y) \ + sup 
It is easy to show 

and 



aij{y) - a. 



\y-m 



Then we have 



Next we consider the C^ "' norms of Dm- We need to recall the proofs of the 
results of [[TOl and [[TTI more carefully. In the case when we consider the L°°- 
norm of Vm for a solution u to the equation (14.61) . the influence of the C^ " norms 
of subdomains appears only in the following constant C in (14.71) : We estimate 
0(|x'|^+°) in the equation (49) in UJj p. 118], i.e. 

fU^') = fUO') + VfmiO')^' + 0(|a;T+") (49) 

as 

|0(|xf +°)| < (4.7) 

(See also flO', Lemma 4.3]). Here C^'" functions /,„ are defined in the cube 
(— 1, 1)", and the graphs of /,„ describe dDm- Now we remark that the constant 
C in (14.71) depends only on the C^'" seminorms of /„. We consider the variable 
change py = x. Then the graph x„ = fm.{x') is changed to ?/„ = fm{y'), where 
fm{y') ■■= p'^fmipy'), and we have 

[/m]ci>"{(-l,l)") < [/m]ci.«((-l/p,l/p)") 

= P"[/m.]ci'"{(-l,l)") ^ [/m]ci'"{(-l,l)")- 

Hence, even when we consider the variable change py = x, we can take the 
constant C in (|4.7I) independent of p. 

Considering the circumstances mentioned above, we can take > inde- 
pendent of p. 
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Now we state the estimate of V^r(x, t; y, s). 
Theorem 4.3. We have 

I V.r(x, t; y, s) I < _ ^^,^,)/, exp (-^f^) (4-8) 

for a.e. x, ?/ G M" and t > s with \x — y\'^ + t — s < 16, where C, c > depend 
only on n, L, jj,, a, A, A, ||oij||c"'("D;;r) seminorms of Dm for some a' 

with (fTTOl) . 

We prove Theorem 14.31 in the same way as the proof of [3 , Proposition 3.6]. 
We first show the following lemmas. 

Lemma 4.4. Le? p := (|xo - + to - r)i/V4- Then 

|r(x,t;^,r)| dxdt<- -—exp' 



for to > T, where C'^, c{ > depend only on n, A, A. 
Proof. By (|4.2I) . it is enough to obtain the estimate 

Jo := / / 77 rrexp ] X[T,oo){t) dx dt 

We consider the following three cases: 

(i) to - < r < to, (ii) to - 2p2 < r < to - p^ (iii) r < to - 2p\ 



Now we consider the case (i). Then we have (vl5 — l)p< |a; — ,^|for any 
X E Bp{xo), because |xo — ^| > a/Ts p. Hence we have 

/o < r / TT^e^P (-7^) dxdt=\B^{0)\p^^ T \,{s)ds, 

where ipi{s) := s""exp(— cip^/s) and ci := 2(a/T5 — l)^c*. If < to — r < 
cip^/n, then we have 

\i{s)ds< £ \,{to-r)ds = {to-rr'^+'exp(^-^^,^ 
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because <^i(s) < (pi{to — r) holds for any s e [0, to — t]. On the other hand, if 
cip^/n < to — r < p^, then we have 

£ ^ v5i(s) ds < £ ^ (^^^ ds = (^^^ {to - T)p-'- exp(-n) 



n\ u _M-n / Cip^ 



< - (io-T)^-"exp 

\ClJ \ to-T 

where we used the properties that 

I for any < s < — r; 

\ n J 

rp' 

n > , and >to — r. 

to-T 

Summing up, we have 

lo < maxjl, [jX\ |Bi(0)|p"(to -r)i-"exp 



Let us consider the case (ii). Then we have (v 14 — l)p < |a; — ^| for all 
X e Bp{xo), because \xo — ^| > V^ip- Hence we have 

/o < /" f j2 — exp [ — _ \ dx dt 

Jtn-o^ JB, 



'to-p^ J Bp{xo) 



{t-rY \ t-T 

= i5i(o)ip'^ r ^2{s)ds, 



where (p2{s) :— s~"exp(— C2p^/s) and C2 := 2(-\/l4 — l)^c*. In a similary way 
as the case (i), if p^ <to —r < C2p^/n, then we have 



/ 

Jta- 



ip2{s) ds < (p2{to -T)ds = p^{to - tY'"' exp I 

Jto-p^-r \ to-T J 

exp (-^). 

because (/92(s) < (p{to—T) for any s G [to—p'^—r, to—r], and we have < to—r. 
On the other hand, if C2p^/n <to — r < 2p^, then we have 

^2{s) ds < ^2 i"-^) ds = (-) " exp(-n) 



<2"-(^y(to-r)-"expf--^ 

C2 / V ^0 - T 
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where we used the properties that 

<P2{s) < ip2 {~~^ any to - - 1" < s < to - 



n>- ,andp > — - — . 

to - r 2 



Summing up, we have 



Jo < |i?i(0)| max |l, 2"-i (^|^ | p"(to - r)^-" exp 



tn — r 



Finally we consider the case (iii). We first remark that 

/•*" I ^—{to- p"^ -T^^^ ifn>2, 

•^^^-p^ i log2 if n = 1, 

because to — r < 2(to — — r). In particular, we have 

r it - r)-" dt < (to - p' - r)-"+i < 2"-i(to - r) 



-n+1 



Hence we have 



/o<i5i(o)ip" {t-Tr-dt<2^-'\B,mp-{to-T)-^+' 

< 2"-^exp(8)|Si(0)|p"(to - r)-"+^xp (-^^^) , 

because |xo - ^\'^/{to - r) < (4p)V2p^ = 8. 

Therefore we have the estimate (14.91 ) in every case. □ 

Now we prove Theorem 14.31 

of Theorem^ Letso,^ G M" and to > r. Letp := (|xo-^P+to-r)i/V4 < 1. 
Then, by Corollary 14.11 we have 

l|Vxr(-, ■;^,^)l|L°°(i?p/2(2^o)x(io-(p/2)2,to)) 

a 

||r(-, ■■,^,T)\\L2(Bpixo)xito-p2,to)), 



pn/2+2 
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because we have 

dV 



i,i=l * ^ 



aij{x)^{x,t;C.,T) ) = Oinfip(xo) x (to -p^^o]• 



By this estimate and Lemma 1441 we have 



II Va;r(-, ■; ^, r) ||Loo(B^/2(^„)x(to_(p/2)2,to]) 
— _9 /J. _\(-«-lW9 ^^P 



p2 (t^_^)(n-l)/2--^ t^_^ 

^ lege: ^ <|xo-eP 

< ^ — exp ' 



(t^ _ r)(n+l)/2 ^ _ ^ 

because we have < 16(to — t)^^- Hence the proof is completed. □ 



Appendix 

In Appendix, we show the estimates (12.11) and (|2.2I) in Lemma [2?T] for the sake of 
completeness. To begin with, we give some embedding lemma which is necessary 
to show the estimates (12.11) and (|2.2I) . First, the following Gagliardo-Nirenberg's 
inequality is well-known (see |I51 p. 24, Theorem 9.3], for example). 

Lemma A.l (Gagliardo-Nirenberg's inequality). Let r, s be any numbers satisfy- 
ing 1 < r, s < oo, and let j, k be any integers satisfying < j < k. If u is any 
function in W,^{W) n then 

||^^m||l9(R") < '^ill^''M||^,(]g„)||M||J^";Jg„), (A.l) 

where 

1 i fl k\ l--f ^^^^ 
- = - + 7 + (A.2) 

q n \s n J r 

for all 7 in the interval 

where a positive constant Ci depends only on n,k,j,r,s,'y, with the following 
exception: If k — j — n/s is a nonnegative integer, then (lA.ll ) holds only for 

j/k < 7 < L 
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Then, as an application of Lemma lA.U we have the following embedding 
lemma. 

Lemma A.2 (embedding lemma). Let Hq{D) be the usual L'^-Sobolev space 
with supports in 13 and v G L"^ {0,T; L^{D)) f] L^ {0,T; H^{D)). Then v e 
^2(n.+2)/nj^Q^ /?o/J5. Moveover, we have the estimate 



|f ||i2(n + 2)/n(Q) < Cl||f 



|2/{n + 2) 

Il°°(0,T;L2(D)) 



\n/(n+2) 



< Ci {\\v\\l^(q^T;LHD)) + \\^v\\l2(Q)) , (A.3) 

where a positive constant Ci depends only on n, and we denote Q := D x (0, T]. 

Proof. We apply Lemma [ATI with q = 2(n + 2) /n, r = 2, s = 2, k = 1 and 
j = 0. Then the equation (IA.2I) yields 'j = n/(n + 2). Hence we have 



Therefore we have 

II ||2(n+2)/n 



< 



2(n+2)/n 

i2(n+2)/n(-£,) "f" 



2(n+2)/n| 



|2/(n+2)\ 



2(n+2)/n 



\lhq)- 

By this inequality and Young's inequality, we have the estimate (IA.3I) . 



dt 



< (^^V"^-^;/"|| ||4/" IIVT 



□ 



Based on Di Giorgi's famous argument, we start to estimate solutions to the 
parabolic equation (11.81 ). By testing max{u — A;, 0}^^ to (11.81 ) we have the follow- 
ing lemma. 

Lemma A.3. Letp > 2. Let Qp := Bp{xo) x (to — p^, to] C Q and ( G C°° ([to — 
p^, to]; (i?p(xo))) satisfy < C < 1 ^^"^ C(') — P^) = 0. T/zen a solution u 
to the parabolic equation (|1.8I) satisfies 

Wiu - A;)+Clli-(to-p2,to;L2(Bp(xo))) + - ||l2(Q^) 



< Co 



[( 


9C 




dt 



L°°{Qp) 



+ llvcilio.(Q^) 1 \\{u-k) 



+ Fl \Qpn{u{x,t) > k} 



l-2/p 



(A.4) 
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for any A; G M, where := max{t>(x), 0}, 

Fo,r '■ = 



(A.5) 



and 6*2 > depends only onn.K and A. 

Proof. Multiplying (11.81 ) by [u — k)^C,'^ and integrating it over Q'p := Bp{xo) 
{to — p^, t') (also see Remark |Z2l ). we have 



X 



(LUS) = jj (^^{u- k)+^ {u- k)+(^dxdt 



_d_ 



{u — k)+ J {u — k)^C'^ dx dt 



2 //q' 



(9 5 



dt 



n 



(9 d 

ain- — ((m — A;)+C)^ — (("^ ~ ^)+C) dxdt 



dx 



dxi 



dx dt 



1 
2 



'Bp {xo) 
n 



{u — fc)^C^ dx 



t=t' 



(u — k)\C— dx dt 
^ '^^ dt 



^ lid (J 

+ 1 1 , ('v-Q^{{u-k)+0-Q^{{u-k)+Odxdt 



d 



J2 a,j{u-k)[ 



dxj dxi 



dx dt. 
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Hence we have 
1 



{u — k)_^_( dx 



Bp{xo) 



t=t' 



+ E 



d d 



i,j=l 



dxj 



dxi 



u — k)\C^ dx dt + '^^^ II aij{u — k) 



Q'p 



dt 



i,j=l 



dxj dxi 



dx dt 



i=l 



d 



+ 11 f{u- k)+Cdxdt + ^ 1 1 fi — {{u-k)+C)dxdt. (A.6) 



We remark that 



d 



fi^{{u- k)+C'^) dxdt 

Q'p 



'pn{u{x,t)>k} 



d 

fi^Q^i('^ - ^)+C) dxdt 



+ 



'pn{u{x,t)>k} 



dc 

fi{u - k)+(—dxdt 

OXi 



Q' r\{u{x,t)>k} 



_d_ 

dxi 



((n-fc)+C) 



dx dt 



1 

+ — 



^1 J JQ'pn{u{x,t)>k} 



IMl^dxdt 



+ 



IfiCl'^dx dt + 

/pn{u{x,t)>k} 

Hence, by (11.41) and (|A.6I) . we have 



'^n{u{x,t)>k} 



{u - kf 



dC 



dxi 



dx dt. 



-[ {u-k)\Cdx \V{{u-k)+C)\^ dxdt 

'^JBpixo) t=t' J J Q'p 



< 



{u — k)X dx 



Bp{xo) 



t=t' 



~i J J O' "^Xo 
l,]=l ^p ■' 
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k)\C^ d,^ + II aijiu — k) 



dt 



i,j=l " " 



' dxj dxi 



dx dt 



+ I l^f(u-k)+Cdxdt + Y^ I I^J^Q^,{iu-k)+e)dxdt 



i=l 



d 



< 



dt 



(M-A;)+da;dt + A||VC||ioc(Q^) // {u-k)Xdxdt 



L^{Qp)J JQ'p 

+ f{u-k)+Cdxdt + ei 



+ I - + 1 

.^1 



'n{u{x,t)>k} 



|V((m - A;)+C)| dxdt 
J2\M'dxdt 



1=1 



+ ^l|vcilioo(Q^) // iu-k)idxdt 



that is 

1 
2 



(u-fc)+C^ dx 



Bp{xo) 

<{A + n) 
+ 



t=t' 



+ // \V{{u-k)+C)\^dxdt 





dC 




dt 



L°^{Qp) 



+ \\^C\\U[Qp) // {u-k)Xdxdt 



1 



y^\f,\'^dxdt+ // f{u-k)+cdxdt. 

Q'pn{u{x,t)>k} .^^ 



Taking the supremum of the inequality over (to — p^, to] with respect to t', we have 

1, 



max 



< (A + n 



2nv«- ^)+Clli-(to-p2,to;L2(Sp(xo)))' A;)+C) 

L°°{Qp) 



i2 



dt 



+ mw^Qp) 



(u — k)'^ dx dt 



- + 1 



Qpn{u{x,t)>k} 



i=l 



J2\fi\^dxdt+ f{u-k)+C^dxdt. 



Qp 



(A.7) 

Now we estimate the last two terms in the right-hand side of (IA.7I ). First we obtain 



pr\{u{x,t)>k} 



\fi\^ dxdt < \Qp n {u{x,t) > k} 
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by Holder's inequality. Now we estimate / f{u — k)^('^ dx dt. We first recall 

fc) + C||L2(n+2)/"(Q,) 

< Ci (\\{U - k) + C\\L^{to-p^,to;LHB,{xo))) + ||V((m - ^) + C) |L2(q^) 



by Lemma IA.21 where Ci > depends only on n. Then, by this inequality, 
Holder's inequality and Young's inequality, we have 



f{u — A;)+C^ dx dt 



L2(n + 2)/n(Q^ 



1 



< e2\\iu - A;) + CllL2(n + 2)/n(Qp) + — ll/CII L2(n + 2)/(n + 4)(Qpn|„(^_t)>;,|) 



<2e2Clmax\\\{u-k)+C\\U^to-p^to-,L-{B,{.o))V || V((m - A;)+C) 



|2 

li^2(Qp) 



+ -||/ll%(™+2) iQpH {u(a;,t) > A;} 

£2 i^^+^(Qp) 



1-2/p 



(A.9) 



because 2(n + 2)/(n + 4) < p(n + 2)/(ra + 2 + By (lATTl) . (ICT) and (ICT) . we 
obtain the estimate (|A.4I) . □ 

By the same argument, we obtain the following lemma for t>_ (x) := max{— t>(a;) 



Lemma IA.3f . Under the same assumption as in Lemma \A.3\ a solution u to the 
parabolic equation (11.81) satisfies 



II (m - A;)_CllL°°(to-p2,to;i2(Bp{a;o))) 



|2 



< Co 



[( 






dt 



+ llvciiio.(Q,) ii(«-fc)-iii2(Q^) 



+ ^o,p|Qpn{u(x,t) < k} 



l-2/p 



(A. 10) 



for any A; G M, where we define F^^p as (|A.5I) . an J C2 > depends only on n,A 
and A. 



The estimate (|2.1I) easily follows from Lemmas I A . 3 1 and [A . 31 . Our next task 
is to prove the estimate (12.21 ). We start by giving a technical lemma which will be 
used later. 
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Lemma A.4. Let C > 0, 6 > 1 and e>O.Ifa sequence {ym}m=o s<^tisfies 

yo < Oo := C-'H-'/'" andO< y„,+, < Ch^y];^^ (A.ll) 

then 



lim ym = Q 

m— >oo 

holds. 

Proof. We show 

Z/m<^, m = 0,l,2,... (A.12) 

by inductive method, where we will determine r > 1 later. By assumption, (IA.12I) 
with m = holds. Hence we now assume (|A.12I) holds, and show (IA.12I) for 
m + 1. By the assumption (lA.l II) and the induction hypothesis, we have 



e 



Now we take r = 6^/^. Then we have 



which is (IA.12I ) for m + 1. □ 

Now we are now ready to show the estimate (12.21 ). The estimate easily follows 
if we have the following lemma. 

Lemma A.5. Let p > n + 2. Then a solution u to (11.81) satisfies the estimate 

\\u\\l^{Qp) < Cp (||m||l2(Q2p) + -^0,2p) , 

where we define Fo,2p by (IA.5I) . and Cp > depends only on n, X, A, p and p. 

Proof. First of all a letter C denotes a general constant depending only on A, A 
andp. Now, let := (l+2-™)p and = A;(2-2-'^) for m = 0,1,2,. . ., where 
we will determine k > later. For m = 0,1,2,..., we take cut-off functions 
Cm e C^iQp J which satisfy 



< Cm < 1 in Q 



Pm 1 
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Cn 



1 in Qp^+i, 



dt 



+ live 



< 



c 



L°°iQp,n) 



m|lL°°(Qp„) ^ / _ \2- 
\Pm Pm+1 ) 



We remark that C™ = on 5p,„(.To) x {to - p^} U dBp^{xo) x (to - p^^o) in 
particular. By Lemmas [A .21 and lA3l we have 

II (m — A;m+l)+Cm||^2(n + 2)/n(Q^^) 

< C(||(m - km+l) + Cm\\loo(^t^_p2^^t^.L2iB,„(xo))) 



< c 



+ 


V{{u- 








[( 


dt 





l|VCm||ioo(Q^^) ) Uu - krn+l)+\\l2(^Qp^) 



+ FL\Qp^r]{uix,t) > kr^+i} 



1-2/p 



< c 



2^"^ 1 2/ 



(A.13) 



where Am{l) ■= Qp^ n {u{x, t) > 1} for / G M. Now we take A; > as 

> pi-{n+2)/p^^^^^_ (A. 14) 

Then we have 

II (m — /i;m+l)+Cm||^2(n + 2)/n(Q^^) 

< c 



2 /c 1 2/ 

— ||(n- A;™+i)+||^2(Q^^^ + ^2(i-(n+2)/p) 

by the estimate (IA.13I) . By defining ip^ '■= II ~ ^m)+|li2(Q^ ), we have 

V^m+l = Wiu — km+l) + Cm\\L2{Qp^^_^) < IK"" — ) + Cm || L2 (Qp„J 

^ |^m(^m+l)| ''IK'^ ^m+l) + Cm||2,2(n+2)/ri('Q^^') 
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22m 



X 



-\\{U - k^+l) + \\L^Q,^) + ^2(l-(n+2)/p) \^rn{km+l] 



+ p2(l-(n+2)/p) \^m{k.m+l)\ 



where we used Holder's inequality and the estimate 

II (m - ^m+l)+|li2(Q^^) < II (m - A;m) + ||^2(Q^^) = <-Pm- 

On the other hand, we have 



V'm = ||(« - fcm)+|lL2(Q,„) > // {u-km)+dxdt 

(fcm+l) 

> I j (km+l - km)X dx dt = 7^^^\A.m{kr, 



that is, 



22m+2 



A;2 



By (|A.15I) and (|A.16I) . we have 



n 4 1+^— o('l ri+2\ ( 4 4 



4 4\ 1 + 



n + 2 p 



We now take k as 



Then we have 



V^m < 



^ > ( 77:; — T / / li^ dxdt\ 
\\Q2p\jJQ2p ] 

// dx dt < // u^dxdt< \Q2p\k 



that is. 



2/p 



< \Q2p?'''k"^. 
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By this inequality and (IA.17I) . we have 



-2 1-^ 



p )k 



,n+2 pj 



< 



> n+2 p> 



n+2 p 



Now we denote := k '^\Q2p\ ^fm- Then by (IA.19I ). we have 
which is the second condition of (|A.11I) with 



(A.19) 



(A.20) 



C = C,h 



i2 H 



^+2 / and e 



n + 2 p 



Then lim^^oo 2/m = if 



yo 



(A.21) 



(A.22) 



by Lemma lA4l where b and £ are defined by (IA.21I ) and C is the constant C in 
(IA.20I) . We remark that the condition (|A.22I ) is equivalent to 



Now we take k as 



{u - k)+\\l2(^Q^^) < eok^\Q2p\ 

1 „ „o 



k' > 



(A.23) 
(A.24) 



Oo\Q2p 

Then the condition (IA.23I ). i.e. the condition (IA.22I) is satisfied. 

Summing up, if we take k such that the conditions (IA.14I ). (IA.18I ) and (IA.24I) 
are satisfied, then we have limm-s>oo Z/m = 0. On the other hand, since 



1 



2p\ 



2p\ 



(w-2A;)+lli2(Q^) 



as m — > oo. 
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Then we have \\{u — 2A;)+||^2(Qp) = 0, that is, 

u<2k a.e. in Qp. (A.25) 

Now we take k as 

1 

a/^o|Q2pI 



which satisfies the conditions (IA.14|) . (IA.18|) and (|A.24I) . Hence we have (|A.25|) . 
which is 

supw < Cp (||u||l2(Q2^) + Fo,2p) . 
Qp 

Replacing Lemma IA.3I by Lemma |A.3( and doing the same argument, we can 
obtain 

-u < Cp {\\u\\lhq^^) + Fo,2p) in Qp 
and thus the proof has been completed. □ 
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